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ABSTRACT 

The inevitable spread in properties of the toroidal obscuration of active galactic nuclei (AGNs) invalidates 
the widespread notion that type 1 and 2 AGNs are intrinsically the same objects, drawn randomly from the 
distribution of torus covering factors. Instead, AGNs are drawn preferentially from this distribution; type 2 are 
more likely drawn from the distribution higher end, type 1 from its lower end. Type 2 AGNs have a higher 
IR luminosity, lower narrow-line luminosity and a higher fraction of Compton thick X-ray obscuration than 
type 1 . Meaningful studies of unification statistics cannot be conducted without first determining the intrinsic 
distribution function of torus covering factors. 

Subject headings: infrared: galaxies - galaxies: active - galaxies: nuclei - galaxies: Seyfert - quasars: general 



1. INTRODUCTION 

The basic premise of the unification scheme is that all 
AGNs are fundamentally the same: accreting supermassive 
black holes. The central engine is surrounded by a dusty 
toroidal structure so that the observed diversity simply reflects 
different viewing angles of an axisymmetric geometry (An- 
tonucci 1993; Urry & Padovani 1995). The classification of 
AGNs into types 1 and 2 is based on the extent to which the 
nuclear region is visible (figure 1). Directions with clear sight 
of the central engine and the broad-line region (BLR) yield 
type 1 sources. Those blocked by the torus result in type 2 
objects, where the existence of the hidden BLR is revealed 
only in polarized light. From basic considerations, Krolik & 
Begelman (1988) concluded that the torus likely consists of a 
large number of individually very optically thick dusty clouds. 
Indeed, VLTI interferometry of the Circinus AGN provide 
strong evidence for a clumpy or filamentary dust structure 
(Tristram et al. 2007). 

In its most extreme form, dubbed the straw person model 
(SPM) by Antonucci (1993), unification posits the view- 
ing angle as the sole factor in determining AGN classifica- 
tion; that is, (1) the classification of any AGN is determined 
uniquely by its viewing angle and (2) the torus is identical 
for all AGNs of the same luminosity. 1 dumpiness invali- 
dates the first assumption because it turns the difference be- 
tween types 1 and 2 into the probability for direct view of the 
AGN (Elitzur 2007, 2008; Nenkova et al. 2008b). Since this 
probability is always finite, type 1 sources can be detected 
from what are typically considered type 2 orientations, even 
through the torus equatorial plane. Conversely, if a cloud hap- 
pens to obscure the AGN from an observer, that object would 
be classified as type 2 irrespective of the viewing angle (see 
figure lc). In cases of such single cloud obscuration, on oc- 
casion the cloud may move out of the line of sight, creating a 
clear path to the nucleus and a transition to type 1 spectrum, 
as observed in a number of sources (see Aretxaga et al. 1999 
and references therein). 

The second assumption underlying SPM, about the torus 
sameness, obviously does not hold. All AGNs cannot be ex- 
pected to have the exact same torus; there must be a spread in 

1 Increasing luminosity tends to decrease the covering factor as first noted 
by Lawrence (1991) in the "receding torus" model. For subsequent observa- 
tions see Maiolino et al. (2007) and references therein. 



torus properties, even among AGNs with the same luminosity. 
Here I discuss some immediate, fundamental implications of 
this inevitable spread. 

2. REALISTIC UNIFICATION 

In panel (a) of figure 1, observer 1 will see the AGN as type 
1, observer 2 as type 2. In panel (b) the AGN orientation is 
the same as in panel (a), only its torus has a smaller covering 
factor. Now both observers see a type 1 object even though 
their viewing angles have not changed. Evidently, the torus 
covering factor Cj is as central to AGN classification as is the 
viewing angle because an AGN with a larger covering factor 
has a higher probability to be viewed as type 2 by a random 
observer. This is obvious also for the more realistic clumpy 
torus, shown in panel (c) of figure 1. Therefore, in a sample 
of AGNs with a distribution of covering factors, those with a 
larger Ct will have a higher probability to be viewed as type-2 
by a random observer, implying that 

AGNs are drawn preferentially from the distri- 
bution of covering factors; type-1 are more likely 
drawn from the distribution lower end, type-2 
from its higher end. 

Contrary to the widespread notion that AGNs of types 1 and 
2 are intrinsically the same objects, fundamental differences 
between their average properties do exist. 

This is a more realistic formulation of the unification 
scheme than SPM, with profound implications for AGN stud- 
ies. Realistic unification immediately explains the findings by 
Ramos Almeida et al. (2009, 2011), based on SED modeling 
of AGN IR emission, that the dusty tori tend to have larger 
covering factors in Seyfert 2 than in Seyfert 1. It also ex- 
plains the seemingly puzzling results of the recent study by 
Ricci et al. (201 1), who analyzed in detail stacked hard X-ray 
spectra (50-200 keV) of all z < 0.2 Seyfert galaxies detected 
with INTEGRAL, the hard X-ray/soft y-ray mission. In agree- 
ment with the basic tenets of unification, both Seyfert 1 and 
Seyfert 2 were found to have the same average nuclear contin- 
uum emission, with a photon index of T = 1.8. But in apparent 
contradiction with unification, the reflection component was 
significantly stronger for the average spectrum of Compton 
thin Seyfert 2 than for Seyfert 1 . Ricci et al find this discrep- 
ancy to arise from a further sub-division among the Seyfert 2 
AGNs. The "lightly obscured" ones (N H < 10 23 cm" 2 ) have 



2 



Moshe Elitzur 





^ 2 





(a) 



(b) 




* 2 



Figure 1. AGN classification in unified schemes, (a) In a smooth-density torus, everyone located inside the cone opening, such as observer 1, will see a type 1 
source; outside — a type 2. (b) Decreasing the torus covering factor, the source becomes a type 1 AGN for more observers, (c) In a clumpy, soft-edge torus, the 
probability for a direct view of the AGN decreases away from the axis, but is always finite. 

the same reflection component as Seyfert 1, R < 0.4, but those 
that are "mildly obscured" (10 23 cirT 2 < A^h < 10 24 cirT 2 ) 
display a much stronger reflection with R - 2.2^] j. While 
this finding contradicts simplistic forms of unification, it is 
precisely the behavior expected from its realistic formulation: 
Seyfert 1 and lightly obscured Seyfert 2 correspond to dif- 
ferent viewing angles of intrinsically similar AGNs, drawn 
from the low end of the covering factor distribution, thus they 
conform, on average, to simplistic unification. But in mildly 
obscured Seyfert 2 the absorber/reflector covers a larger frac- 
tion of the X-ray source, producing stronger reflection that is 
not seen in the average Seyfert 1 spectrum, where large cover- 
ing factors are rare. The large difference between the average 
reflection spectra of Seyfert 1 and 2 arises from significant 
differences in their torus covering factors. 

3. COVERING FACTORS 

The sometime loosely invoked concept of "torus covering 
factor", Ct, can be rigorously defined as the fraction of the 
sky at the AGN center covered by obscuring material; it is the 
same as the fraction of randomly distributed observers whose 
view to the center is blocked, and thus see a type 2 AGN 
(Nenkova et al. 2008a). Denote by N(i) the overall number 
of clouds encountered, on average, along angle i from the 
axis. Then the probability for direct viewing of the AGN 
from that direction is e~ N ^ and the torus covering factor is 
Cj = 1 - J e~ N ^d cos i. If A^o is the average number of 
clouds along radial equatorial rays then N(i) = No<3}(i), where 
<E>(90°) — 1. The cloud angular distribution function O can 
be conveniently parametrized as Gaussian, fl>(z) = e~( 90 ~')"/°" , 
with cr the distribution angular width (Elitzur et al. 2004; 
Nenkova et al. 2008b). Fitting of IR observations with clumpy 
torus models with Gaussian angular distributions has been re- 
ported by a number of teams, and figure 2 shows the results 
of these modeling efforts in the No-cr plane together with the 
contour plots of Ct- 2 As expected from realistic unification, 
and first noted by Ramos Almeida et al. (2009), type 1 and 
type 2 AGNs preferentially occupy different regions in the 
plane. The few sources with a "wrong" covering factor (large- 
Ct type 1, small-CT type 2) merely reflect the probabilistic na- 
ture of clumpy obscuration. Although this ad-hoc collection 

2 Earlier versions of this figure were presented in Elitzur (2009; acces- 
sible at www.mpe.mpg.de/events/pgnS9/online_proceedings.html) 
and Ramos Almeida et al. (201 1). The contour plots in both of these earlier 
figures are afflicted by the computer bug reported in Nenkova et al. (2010). 




Figure 2. Clumpy torus covering factors. Contour plots are for a toroidal 
Gaussian distribution, where the number of clouds along viewing angle i from 

the axis is A'oe~ (90 ~ ! '"'' " , with No and cr free parameters. Each contour is the 
locus of Nq-<t combinations that produce the labeled covering factor. The 
data points are from clumpy torus modeling of IR observations of AGNs 
reported in Mor et al. (2009), Nikutta et al. (2009), Ramos Almeida et al. 
(201 1), Alonso-Herrero et al. (201 1) and Deo et al. (201 1). 

of AGNs, which were selected by different, unrelated crite- 
ria, does not constitute a complete sample (only the Mor et al. 
2009 analysis of PG quasars involved a complete sample), it 
does illustrate the point. 

Since the covering factor measures the fraction of AGN lu- 
minosity captured by the torus and converted to infrared, the 
AGN IR luminosity is CjL, where L is its bolometric lumi- 
nosity. Therefore type 2 AGNs have intrinsically higher IR 
luminosities than type 1. Contrary to earlier expectations of 
strong anisotropy at A < 8/mi, Spitzer observations at this 
wavelength regime show a great similarity between the IR 
fluxes of AGNs 1 and 2 when normalized with either their 
X-ray fluxes (Lutz et al. 2004; Horst et al. 2006) or optically 
thin radio emission (Buchanan et al. 2006). Part of this puzzle 
was solved by clumpy torus calculations, which show much 
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less anisotropy in IR emission than the earlier smooth-density 
models (Nenkova et al. 2008b). Realistic unification explains 
away the remainder. 

While infrared arises from reprocessing of the AGN radia- 
tion captured by the torus, narrow line emission is generated 
by the radiation that has escaped the torus. The narrow line 
luminosity is proportional to (1 - Ci)L. At the same bolomet- 
ric luminosity, type 1 AGNs can be expected to have a higher 
narrow line luminosity than type 2. 

4. UNIFICATION STATISTICS 

Implicitly or explicitly, all studies of AGN statistics assume 
that type 1 and type 2 are intrinsically the same objects, drawn 
randomly from the distribution of torus covering factors. Re- 
alistic unification invalidates this assumption, implying that 
covering factors deduced from analysis of type 2 sources do 
not necessarily apply to type 1 and vice versa. This may 
contribute to discrepancies among studies of AGN statistics. 
Thus the fraction of Compton thick AGNs (X-ray obscuring 
column A^h <: 10 24 cirT 2 ) was determined to be as high as 50% 
in a pre-selected sample of Seyfert 2 (Risaliti et al. 1999; 
Guainazzi et al. 2005) but only ~ 20% in complete X-ray sam- 
ples without spectroscopic pre-selection (Malizia et al. 2009; 
Burlonetal. 2011). 

Schmitt et al. (2001) find that the type 2 fraction among 
Seyfert galaxies is ~70%, while Hao et al. (2005) find it to 
be only ~50%. This discrepancy can be attributed, at least in 
part, to different selection methods. Both studies based their 
selections on what they considered reliable isotropic proper- 
ties: Schmitt et al used 60 //m IR emission, Hao et al the 
narrow-line luminosity of [O in] /15007. However, as noted 
above, thanks to the larger covering factor of their torus, 
Seyfert 2 convert a larger fraction of their luminosity to IR, 
therefore IR-selection is biased in their favor. The opposite 
afflicts the line selection criterion of Hao et al, which in- 
troduces preference for smaller covering factors and type 1 
AGNs. Even though the emission itself might be isotropic, 
in both cases it involves reprocessing of different fractions of 
the intrinsic UV/optical continuum, and these fractions intro- 
duced opposite biases in the two studies. 

All previous findings involving unification statistics, in- 
cluding the synthesis of the cosmic X-ray background (e.g., 
Gilli et al. 2007; Treister et al. 2009), therefore need a critical 
reexamination and revision. One cannot draw statistical in- 
ferences from AGN populations without folding in the intrin- 
sic distribution of torus covering factors, which is unknown. 
A reliable determination of this distribution function requires 
IR modeling of an unbiased, complete sample of AGNs. The 
only practical method to select AGNs with minimal classifica- 
tion bias is through hard X-ray surveys; although such selec- 
tion still misses the most heavily obscured (Nh > 10 25 crrT 2 ) 
Compton thick AGNs (Burlon et al. 2011; de Rosa et al. 
2012), it is the least biased method of AGN selection. From 
the catalogue of INTEGRAL observations in the 20^-0 keV 
band, Malizia et al. (2009) extracted a complete sample of 88 
AGNs, including type classification. Using 0.3-195 keV data 
from Swift-BAI observations and the XMM-Newton archive, 
Burlon et al. (2011) have compiled another complete sample 
that contains 199 type-classified AGNs. Identifying Spitzer 
counterparts to AGNs in these samples and fitting their IR 
observations with clumpy torus models, as was done for the 
sources shown in figure 2, is the only feasible approach to 
determining the intrinsic distribution function of torus cov- 
ering factors. Searching through the Spitzer archives shows 



that the Malizia et al sample has 54 matches (L. Bassani and 
A. Malizia, private communication); a similar search for the 
Swift-BAI sample can be expected to add significantly to the 
number of blindly selected AGNs with Spitzer counterparts. 
In a future publication we will report separately the results of 
clumpy torus modeling of the IR emission for each source in 
this expanded list to determine its torus covering factor from 
the fitted values of A^o and cr, and construct from the entire 
sample the probability distribution function of Cj- 

5. DISCUSSION 

Introducing SPM, Antonucci (1993) suggested that "The 
true situation may be in some sense half way between the 
SPM and the hypothesis that orientation doesn't affect clas- 
sification at all," ascribing classification to a single property 
of the system — orientation. However, quite apart from the 
indeterminism introduced by dumpiness, AGN classification 
involves a fundamental plane spanned by two independent 
axes — orientation and covering factor. Simply put, the true 
situation is not somewhere between SPM and no orientation 
effects; it involves an additional, independent variable. 

Because the covering factor is an intrinsic, observer- 
independent property, it makes what has always been a dif- 
ficult problem even more difficult. AGN selection by an 
isotropic emission indicator is insufficient for reliable statis- 
tics because most measured fluxes involve reprocessing of the 
intrinsic UV/optical radiation and the fraction captured for re- 
processing may differ. The IR luminosity is C T L while the 
narrow line luminosity is proportional to (1 - Ci)L, introduc- 
ing biases in favor of, respectively, types 2 and 1 . Past studies 
of unification statistics cannot be fully trusted before they are 
repeated taking into account the, as yet unknown, distribution 
of covering factors. 

I thank Loredana Bassani and Angela Malizia for help with 
their AGN sample. Support by NASA and NSF is gratefully 
acknowledged. 
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